翻訳と辞書
Words near each other
・ Split-Hopkinson pressure bar
・ Split-horizon DNS
・ Split-intein circular ligation of peptides and proteins
・ Split-Karlovac Fault
・ Split-level hogfish
・ Split-level home
・ Split-octonion
・ Split-off point
・ Split-phase electric power
・ Split-pi topology
・ Split-quaternion
・ Split-radix FFT algorithm
・ Split-rail fence
・ Split-ring resonator
・ Split-single
Split-step method
・ Split-T
・ Split-ticket voting
・ Split/EP
・ Split/Second
・ Splitboard
・ Splitfin flashlightfish
・ SplitFire spark plug
・ Splithand Township, Itasca County, Minnesota
・ Splitlog Church
・ Splitrmx12
・ Splitrock (Austin, Texas)
・ Splitrock Reservoir
・ Splitrock, Minnesota
・ Splitska


Dictionary Lists
翻訳と辞書 辞書検索 [ 開発暫定版 ]
スポンサード リンク

Split-step method : ウィキペディア英語版
Split-step method
In numerical analysis, the split-step (Fourier) method is a pseudo-spectral numerical method used to solve nonlinear partial differential equations like the nonlinear Schrödinger equation. The name arises for two reasons. First, the method relies on computing the solution in small steps, and treating the linear and the nonlinear steps separately (see below). Second, it is necessary to Fourier transform back and forth because the linear step is made in the frequency domain while the nonlinear step is made in the time domain.
An example of usage of this method is in the field of light pulse propagation in optical fibers, where the interaction of linear and nonlinear mechanisms makes it difficult to find general analytical solutions. However, the split-step method provides a numerical solution to the problem.
==Description of the method==
Consider, for example, the nonlinear Schrödinger equation
: = - + i \gamma | A |^2 A = (D + \hat N )A,
where A(t,z) describes the pulse envelope in time t at the spatial position z. The equation can be split into a linear part,
: = - = \hat D A,
and a nonlinear part,
: = i \gamma | A |^2 A = \hat N A.
Both the linear and the nonlinear parts have analytical solutions, but the nonlinear Schrödinger equation containing both parts does not have a general analytical solution.
However, if only a 'small' step h is taken along z, then the two parts can be treated separately with only a 'small' numerical error. One can therefore first take a small nonlinear step,
:A_N(t, z+h) = \exp\left(\gamma |A|^2 h \right ) A(t, z),
using the analytical solution.
The dispersion step has an analytical solution in the frequency domain, so it is first necessary to Fourier transform A_N using
:\tilde A_N(\omega, z) = \int_^\infty A_N(t,z) \exp() dt ,
where \omega_0 is the center frequency of the pulse.
It can be shown that using the above definition of the Fourier transform, the analytical solution to the linear step, commuted with the frequency domain solution for the nonlinear step, is
:\tilde(\omega, z+h) = \exp\left(inverse Fourier transform of \tilde(\omega, z+h) one obtains A\left(t, z+h\right); the pulse has thus been propagated a small step h. By repeating the above N times, the pulse can be propagated over a length of N h.
The above shows how to use the method to propagate a solution forward in space; however, many physics applications, such as studying the evolution of a wave packet describing a particle, require one to propagate the solution forward in time rather than in space. The non-linear Schrödinger equation, when used to govern the time evolution of a wave function, takes the form
:i \hbar = - + \gamma | \psi|^2 \psi = (D + \hat N )\psi,
where \psi(x, t) describes the wave function at position x and time t. Note that
:\hat D=- and \hat N =\gamma | \psi|^2 , and that m is the mass of the particle and \hbar is Planck's constant over 2\pi.
The formal solution to this equation is a complex exponential, so we have that
: \psi(x, t)=e^\psi(x, 0).
Since \hat and \hat are operators, they do not in general commute. However, the Baker-Hausdorff formula can be applied to show that the error from treating them as if they do will be of order dt^2 if we are taking a small but finite time step dt. We therefore can write
: \psi(x, t+dt) \approx e^e^\psi(x, t).
The part of this equation involving \hat N can be computed directly using the wave function at time t , but to compute the exponential involving \hat D we use the fact that in frequency space, the partial derivative operator can be converted into a number by substituting ik for \part \over \part x , where k is the frequency (or more properly, wave number, as we are dealing with a spatial variable and thus transforming to a space of spatial frequencies—i.e. wave numbers) associated with the Fourier transform of whatever is being operated on. Thus, we take the Fourier transform of
:e^\psi(x, t),
recover the associated wave number, compute the quantity
: e^,
and use it to find the product of the complex exponentials involving \hat N and \hat D in frequency space as below:
: e^F(Fourier transforms of this algorithm can be computed relatively fast using the ''fast Fourier transform (FFT)''. The split-step Fourier method can therefore be much faster than typical finite difference methods.

抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)
ウィキペディアで「Split-step method」の詳細全文を読む



スポンサード リンク
翻訳と辞書 : 翻訳のためのインターネットリソース

Copyright(C) kotoba.ne.jp 1997-2016. All Rights Reserved.